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Convective Core Size Determination Using Asteroseismic Data 
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As it is long known, the presence of a convective region creates a discontinuity in the chemical profile of a star, which in 
turn translates into a sharp variation of the adiabatic sound speed. This variation produces an oscillatory behavior of the 
pulsation frequencies, to which low degree p-modes are sensitive. We investigate the possibility of detecting the signature 
related to the presence of a convective core in the frequency spectrum of low-mass stars by means of suitable frequency 
combinations (such as separations and ratios). 
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1 Introduction 

Asterosesimology has shown to be the most powerful! tool 
to study the inner properties of stars through the observation 
of stellar oscillations. Many efforts have been conducted 
to understand the interior physics of stars by understand- 
ing the influence that different physical phenomena have on 
the frequency spectra. In particular, it is known that sharp 
variations on the adiabatic sound speed produce an oscil- 
latory signal in the p-mode spectra whose period is related 
to the position of such a discontinuity (e.g.'V orontsov 19881 
lAudard & Provost I994| l. Bearing this in mind, in this study 
we apply different asteroseismic tools to a computational 
model of a star containing a convective core and try to ex- 
tract its position by means of the expected signatures on the 
frequency spectra. 



^ 2 Signature of the convective core 



During its main sequence evolution stars with masses higher 
than '^l.l Mq develop a convective core, whose extension 
is determined by the temperature stratification on the inte- 
rior. Since the timescale for mixing of elements in this con- 
vective region is much shorter than the nuclear timescale, 
the core is believed to be homogeneously mixed in compar- 
ison with the radiative envelope, and a discontinuity in den- 
sity appears in the edge of the fully mixed core. This sharp 
variation in density translates into a discontinuity in the adi- 
abatic sound speed, to which p-modes are sensitive. Theo- 
retically, the imprint left in the frequency spectrum of the 
acoustic modes is the appearance of an oscillatory signal, 
which is related to the location of the region of sharp varia- 
tion within the star In particular, the period of this oscilla- 
tion relates to the exact radial coordinate where the disconti- 



nuity is located, and therefore to the travel time of the propa- 
gating wave throughout the acoustic cavity dMonteiro et al. 1994] l.[ 
The travel time of a propagating wave is represented by the 
acoustic radius or the acoustic depth, which are alternative 
representations of each other and are given by: 



, dr 

/o Cs 



dr 



(1) 



where Cg is the adiabatic sound speed. If the location of the 
density discontinuity is given by, say, ri in radial coordi- 
nates, and that same position is represented in acoustic ra- 
dius and acoustic depth by ti and ri respectively, the period 
of the oscillation induced by this sharp variation is given by 
l/(2ti) and 1/(2ti). 

3 Frequencies, separations, and ratios 

In order to study the inner structure of a star using astero- 
seismology, we need to find the appropriate tool which al- 
lows us to probe the desired region of the star an extract the 
required information. Different combinations of low-degree 
p-modes have been suggested as suitable probes of the phys- 
ical characteristics of a star (e.g. Christensen-Dalsgaard 1984) l.[ 
The most commonly used are the so caUed large and small 
frequency separations, defined as: 



Ai{n) = - 1^71-1,1 



(2) 
(3) 
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where Vn,! is the mode frequency of angular degree I and 
radial order n. However, these combinations are affected by 
the outer layers of the star where turbulence in the near- 
surface is still poorly understood, and also changes in the 
adiabatic sound speed occur due to the presence of the con- 
vective envelope and the helium second ionization zone (e.g.| 
IBallot et al. 20041 l. For the purpose of our study, we need to 
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isolate from surface contamination the signal arising from 
the interior of the star in order to properly quantify the ef- 
fects in the frequency spectra of a discontinuity in the chem- 
ical profile outside of the convective core. Therefore, we 
consider here the smoother 5 points small frequency sepa- 
rations and the ratio of small to large separations, quantities 
which are mainly determined by the inner structure of the 
star ( [Roxburgh & Vorontsov 2003) and are constructed as: 



'4i^„_i4-^6j/„,o-4j^n,i + ^'n+i,o) (4) 



dw = -g(i^«-i,i-4z^„_o + 6:^„,i-4i/„+i,o + J^„+i,i) (5) 



dw{n) 
Ao(n+ !)■ 



(6) 
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The small frequency separations have already been used to 
identify the location of the convective envelope and the he- 
lium second ionization zone (Hell) in the sun by applying 
them to observational data ( [Roxburgh 2009 1, while it has 
been shown that the ratios fairly cancel out the influence 
in the frequencies of the outer layers ( Roxburgh 2005| l. We 
will explore here their potential in a more massive model 
than the sun to extract the signature produced by the sharp 
variation in the adiabatic sound speed outside of the convec- 
tive core. 



4 Model computation 

We have used a 1.5 Mq model at solar metallicity con- 
structed with GARSTEC (Garching Stellar Evolution Code, 
IWeiss & Schlattl 2008l l during its main sequence evolution. 
Its central hydrogen content is Xc'^ 0.25 and nuclear burn- 
ing has produced a discontinuity in the chemical composi- 
tion outside the fully mixed convective core. The adiabatic 
sound speed profile of the model is shown in Fig.lT] where 
the location of the convective core, the convective envelope 
and the Hell ionization zone are also depicted. We calcu- 
lated the frequencies for the model using ADIPLS (Aarhus 
Adiabatic Oscillation Package, Christensen-Dals gaard 2008| [ 
and constructed the small separations and ratios defined in 
Eqs.[4]|5]and|6] Both quantities are plotted in Fig. [2] where 
a clear long period oscillatory component is observed. As 
explained in Sect. [2] the period of this oscillatory imprint 
in the frequency spectra can provide us with information 
about the location of the sharp discontinuity producing it, 
so we have done a non-linear least squares sinusoidal fit to 
the ratios and differences, and calculated the period of the 
obtained sine wave. The result of these fits are also plotted 
in Fig. in where the value obtained for the period of each 
curve is given on the upper-right side of each panel. We find 
that the long period oscillation observed is produced at a po- 
sition of approximately ^187 s for the small frequency sep- 
arations and at ~134 s for the frequency ratios, both values 
being in the vicinity of the convective core location (~135 
s, see Fig.dj). The results are summarized in Tabled] where 



Fig. 1 Adiabatic sound speed profile for the 1.5 M© 
model (solid black line). Also plotted in blue dashed lines 
are the location of the convective core boundary (Ccore), 
the base of the convective envelope (BCE) and the helium 
second ionization zone (Hell). 
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T[s]= 134. 4 
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Fig. 2 Small separations (upper panel) and ratios of small 
to large separations (lower panel). The frequencies marked 
in red have been used for the sine wave fitting (black solid 
line), and the period of the sine wave with its corresponding 
acoustic radius are also shown. The vertical blue dashed line 
corresponds to the acoustic cut-off frequency. The value of 
the period and its corresponding acoustic radius are on the 
upper-right side of each panel. See text for details 



the estimated position of the convective core extracted us- 
ing the seismic tools is shown in mass and radius coordi- 
nates. We see that the ratios provide a more accurate loca- 
tion than the small separations, probably because the small 
separations still retain information from the outer parts of 
the star which is not being completely cancelled by this fre- 
quency combination (as it was already shown for the sun, 



see|Roxburgh2009l. 
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Table 1 Convective core location - in terms of acoustic 
radius (in seconds), fraction radius and fraction mass - from 
the model, and the estimated location from the frequency 
combinations. 





t 


m(r)/R 


r/R 


Model 


135.0 


0.065 


0.058 


ddoi & ddio 


187.0 


0.132 


0.081 


roi & rio 


134.4 


0.065 


0.058 



5 Turning points 

The information about the interior of the star that each acous-l 
tic mode carries depends in the region of the start it probes, 
and therefore in the acoustic cavity where it propagates. The 
radial (l — 0) modes propagate down to the center of the 
star, but the situation for the I = 1 modes is not so simple. 
In a first approximation, and neglecting the perturbation in 
the gravitational potential (Cowling 1941), the propagation 



cavity of the p-modes can be estimated for each frequency 
value and the inner turning point calculated. This is the so 
called characteristic acoustic frequency, given by: 



Cs{rt)L 



rt 



(7) 



where rt is the radius where the turning point is located, 
and L = (Z + 1/2) for I = 1 modes. When the gravitational 
potential is taken into account, the fluctuations on it modify 
the dispersion relation of the incoming wave lJeans 1928l and 
the characteristic frequency is corrected as: 



47rGp(rt), 



(8) 



where p{rt) is the density in rt- The inner turning point 
for each frequency value corresponds to the radius at which 
Eq.|7]or Eq.|8]is satisfied. We calculated this values for the 
I ~ 1 modes in both formulations to check if the modes 
used in the sinusoidal fit shown in Fig. |2] actually penetrate 
the convective core. The obtained characteristic frequencies 
are plotted in Fig. [3] as a function of radius. The green line 
represents the prescription from Eq. |71 while the blue solid 
line shows the result for Eq.|8] It can be seen in Fig.|2]that 
there is a clear change in the behavior of the small sepa- 
rations (and the ratios) around ^ 800 ^Hz. In the frame 
of the Cowling approximation, this cannot be explained: 
only dipolar modes above ^ 1700 /iHz should reach the 
core. However, by considering the fluctuations in the grav- 
itational potential (Eq. |8]l, I = 1 modes above ~ 800 /iHz 
can reach the core - even if it is worth noting that below 
^ 1200 /iHz, they must cross an acoustic wall (see Fig.[3]i. 
Thus this corrected Lamb frequency appears to be a suited 
approximation to know which modes can reach the core or 
not. For a future exploration, we will considered the results 
of ITakata 20061 who has developed a self-consistent theory 
for dipolar modes without the Cowling approximation. 
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0.16 



Fig. 3 Characteristic acoustic frequency for the I = 1 
modes. Blue line shows the calculation made including the 
gravitational potential, while the green line depicts the one 
ignoring it. Black dashed line represents the position of the 
Convective Core, while the solid red line shows the fre- 
quency value of the lowest I = 1 mode used to perform 
the sinusoidal fit described in Sect.|4] See text for details. 

6 Conclusions and future development 

The diagnostic potential of the long known periodic com- 
ponent of frequency perturbations due to a discontinuity or 
steep gradient in the adiabatic sound speed has been used 
to estimate the position within the acoustic cavity of such 
sharp feature. In particular, we have estimated the location 
of the convective core using a sine wave fit to the large os- 
cillation pattern observed in two different frequency com- 
binations, showing that the accuracy in the position deter- 
mination is higher when the ratios are considered than us- 
ing the small separations. We intend to continue this study 
and extend it to other mass ranges and different types of 
discontinuities, such as the ones produced by overshoot- 
ing and semiconvection ( Silva Aguirre et al. 2010[ ). We are 
aware that the potential of this diagnosis depends on the 
ranges of radial orders considered for the sine wave fitting of 
the data, for which a more elaborate analysis of the turning 
points of the dipole modes will be applied (e.gl Takata 20061) .[ 
Finally, we plan to complement this study by also apply- 
ing other asteroseismic tools developed to detect small con- 
vective cores to our models (e.g ICunha & Metca lfe 20071 
IMazumdar et al. 2 006 ), and explore the possibilities of de- 
tecting these mixed zones with the current data being ob- 
tained by the CoRoT and Kepler missions. 
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